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DIRECT EVIDENCE FOR TERMINATION OF OBSCURED STAR FORMATION BY RADIATIVELY DRIVEN 

OUTFLOWS IN REDDENED QSOS 

Duncan Farrah\ Tanya Urrutia^, Mark Lacy^, Andreas Efstathiou^, Jose Afonso^'^, Kristen Coppin^, Patrick 
B. Hall^, Carol Lonsdale^, Tom Jarrett^, Carrie Bridge^°, Colin Borys^, Sara Petty^^ 

ABSTRACT 

We present optical to far-infrared photometry of 31 reddened QSOs that show evidence for radia- 
tively driven outflows originating from AGN in their rest-frame UV spectra. We use these data to 
study the relationships between the AGN-driven outflows, and the AGN and starburst infrared lu- 
minosities. We find that FeLoBAL QSOs are invariably IR-luminous, with IR luminosities exceeding 
lO^^L© in all cases. The AGN supplies 76% of the total IR emission, on average, but with a range 
from 20% to 100%. We find no evidence that the absolute luminosity of obscured star formation is 
affected by the AGN-driven outflows. Conversely, we find an anticorrelation between the strength of 
AGN-driven outflows, as measured from the range of outflow velocities over which absorption exceeds 
a minimal threshold, and the contribution from star formation to the total IR luminosity, with a much 
higher chance of seeing a starburst contribution in excess of 25% in systems with weak outflows than 
in systems with strong outflows. Moreover, we find no convincing evidence that this effect is driven 
by the IR luminosity of the AGN. We conclude that radiatively driven outflows from AGN can have 
a dramatic, negative impact on luminous star formation in their host galaxies. We find that such 
outflows act to curtail star formation such that star formation contributes less than ~ 25% of the 
total IR luminosity. We also propose that the degree to which termination of star formation takes 
place is not deducible from the IR luminosity of the AGN. 

Subject headings: accretion, accretion discs, quasars: absorption lines, galaxies: starburst 



1. INTRODUCTION 

The last decade has seen substantial progress in un- 
derstanding how galaxies assemble their stars and cen- 
tral supermassive black holes (SMBHs). Recent results 
almost all point to the same conclusion; that the bulk 
of galaxy assembly occurred at z > 0.7, and that a 
significant fraction of it took place in obscured 'bursts' 
of intense star formation and SMBH accretion. Indi- 
rect evidence for this comes from, for example, the ste l- 
lar ages of low-redshift galaxies (iHeavens et al.l 120041 ) , 
the ex istence of evolved ehiptical galaxies at high red- 
shifts (jPunloD et all 119961 : lElIis et al.l [19971 : iRakos et all 
120071). and studies of stellar mass assembly over wide 
redshift ranges CG lazebrook et al. 2004; Fontana et al. 
120061 : iPerez-Gonzdlez et al.ll2008l : iMarchesini et al.l 120091 : 
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lllbert et al.ll201Qh . Direct evidence mostly comes from 
extragalactic infrared and millimetre imaging surveys, 
which find that the number density of IR-luminous 
galaxies i ncreases dramatically with increasing red- 
shift (e.R. iLe Floc^h efal] l2005UC hapmai renn 
Babbedge et all 120061: ICoppin e t al. 2006: iShupe et al.l 



20081 : lAuster mann et al.1 120091 : Bales e t all 1201 Of ), and 



that mu ch of the growth period of SMBHs was shrouded 
in dust ([Martmez-Sansigre et al.l 120051 : 1 Alexander et al.l 
120081) . 

The importance of obscured starburst and AGN ac- 
tivity in assembling galaxies suggests that the two 
phenomena may affect each other. A link between 
them is implied by, for example, the tight correla- 
tion between the mass of the SMBH and both stellar 
bulge mass (Magorrian et al. 1998; Ferrarese & Merritt 
1 2006: Tremaine et al.1 120021: Marconi & Hunt 2003; 
IShields et al . 2003; iHaring fc Rix 2 004) and dark mat- 
ter halo mass (jFerraresd l2002[ ). and the coeval 
presence of both starburst and AGN activity in 
many IR-luminous systems at all redshifts (see e.g. 
Blain et al. 2002; Lagache et al. 2 0051 : iLonsdale et al.l 
2006; Hernan-Caballero et al. 2009). 

Recently, interest into the relationship between star- 
burst and AGN activity has been stimulated by what first 
appeared to be irreconcilable differences between obser- 
vational and theoretical results on the assembly history 
of galaxies. Foremost among these were the difficulties 
that models faced in explaining the observed galaxy lumi- 
nosity function (LF) at low and high redshifts simultane- 
ously; if the models were tuned to match the local galaxy 
LF then they underpredicted the number of massive 
galaxies observed at high redshifts, whereas the models 
that matched the number of high redshift galaxies gave 
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The Sample. 
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Note. — Positions, redshifts and SDSS magnitudes are taken from the SDSS DR6. For the SDSS magnitudes we used 
the PSF magnitudes in the AB system. For 2MASS, we used the 'default' magnitudes in the Vega system, taken from the 
public 2MASS All-Sky Point Source Catalogue. WISE magnitudes from the operational database as of June 2011. An ' ' 
indicates that the source is not in the catalogue in that band. 
^ UKIDSS LAS magnitude 



poor fits to the l o cal galaxy LF (iKauffmann et al.l ll999l: 
Cole et alJ I2QQQI : iSomerville et all I2QQ1I : iBens on et alJ 
2QQ3I ). Other problems included: (1) the prediction 



that a higher fraction of rich galaxy cluste rs should har- 
bour cooling flows than is observed (Pete r son et alJl2QQ3l : 
IXu et alJl2QQllPeterson fc Fabianll2QQ6f ). (2) the diflicul- 
ties that early semi- analytic models faced in reproducing 
the large number of IR- luminous g alaxies observed at 
high redshift (e.g. Baugh et al. 2005), and (3) the expec- 
tation that the mass return rate from stars would cause 
central SMBHs to be larger than is observed (Ciotti et J] 
[19911 : ^Ciotti fc Qstriker 2007) . 

Several solutions have been proposed to resolve these 
issues. Among the most promising is an idea termed 
'AGN feedback'. AGN feedback is the exertion of influ- 
ence of an SMBH on kpc to Mpc scales to curtail star 
formation in the host galaxy, and/or further accretion 
onto the SMBH itselQ This can occur in four ways. 
Radiation from regions immediately local to the SMBH 
can (1) heat gas in the ISM so it cannot collapse to form 
stars, and/or (2) radiatively drive gas out of the galaxy 
(or, equivalently, stop it falling in from the IGM), thus 
emptying the galaxy of fuel for further star formation. 
Additionally, matter expelled from regions immediately 

^'^ The idea that SMBHs may exert 'negative' feedback on in- 
flowing gas to regulate their own growth and that of the host galaxy- 
predated its use to reconcile galaxy evolution m odels with obser- 
vations, see e.g. Haiman et al. 199^ ; ISilk ReedflQQS 



local to the SMBH can (3) heat gas in the ISM, and/or 
(4) drive gas out of the galaxy via kinetic pressure. 

Most discussions of AGN feedback condense the four 
feedback mechanisms described above into two simpli- 
fled paradigms; 'quasar' mode feedback and 'radio' mode 
feedback. Quasar mode feedback assumes that radia- 
tion from an accretion disk terminates star formation 
in the host galaxy, usually by coupling some fraction of 
the QSO luminosity to kinetic energy injected into the 
ISM. Quasar mode feedback is thought to be short-term, 
acting for only the ~ 10^ years over which the quasar 
is accreting rapidly. Radio mode feedback on the other 
hand is assumed to occur via a relativistic jet that trans- 
fers momentum to the ISM. Radio mode feedback can 
occur over longer timescales (~ 10^ years) than quasar 
mode as it requires only a low accretion rate to produce 
jets sufliciently powerful to affect the ISM. Nevertheless, 
radio mode and quasar mode feedback are not mutually 
exclusive. 

For individual galaxies, AGN feedback is predicted 
to have a significant impact, both from analytic 
(Fabian _1 999; Wvithe & Loeb 2003 ; iSazonov et al.l 



2005; Fabian et al.l 120061: iPopd 120091 : IPower et al.l 120111: 
Kavirai et al. 112011') and num erical studies (|O mma et al.J 
2004; Ciotti & Ostriker 200 71: IAntonuccio-l5e logu & Silk' 
[ 2008.: .To rtora et al.i |2009l: lOstrike r et al. .2010i : 
IWamer Bicknell l2C)lll : iHambrick et al. 2011). For 
galaxy mergers, quasar mode feedback is predicted 
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TABLE 2 

MIPS FLUXES, ABSORPTION STRENGTH MEASURES (SEE EQUATION [T] & TABLE [S]), AND IR LUMINOSITIES. 
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Infrared Luminosities (log (L©)) 
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Note. — Errors on the absorption strengths from uncertai nties in the fits are 1 00-200 km s ^, but see ^4.21 fc the appendi x. The AS^^^ 
values are an alternative measure of ASq (Table [3]), taken from [Gibson et al.| (120091) (a ' ' indicates that the object is not in the lGibson et all 
(I20Q9D sample). Errors are the 90% confidence intervals, see ^3.21 
^ Total infrared luminosity, integrated over l-lOOO/im in the rest-frame. 
^ Infrared luminosity of the AGN component, integrated over l-lOOO/xm in the rest-frame. 
^ Infrared luminosity of the starburst component, integrated over 1-1000/im in the rest-frame. 



to have a profQund influence (iDi Matteo et all I2QQ51; 
Springel et all I2QQ5I : iDebuhr et all I2Q11I : iSnvder et all 
2Qllf ). Models for the cosmological evolution of 
galaxies and clusters have incorporated one or both 
of these feedback paradigms, often with noticeable 
improvements in reproducing observations. These 
i nclude sem i-analytic models using quasar mode 
Granato et a L 2004; iMenci et all 120081), radio mod e 
Bower et al. 2006, 2008; Gonzalez-Perez et alJ "2009"), 
or both (Somerville et al. 2008), and numer ical models 
with quasar mo de (IBooth fc: Scha ve 200 91; ISales et al.l 
20T0I: iMcCa rthv et all 120111 : iTevssierlFldl [20TTI : 
Chatteriee et al. 201 1|), or both forms of feedback 
(iCroton e t al. 2006; De Lucia et al. "2006"; "Siiacki et alj 
20071: [Kitzbichler & White 2007; Monaco et al] 
20071: Fu chwein et all 120081 : iShort &: ThomasI 120091 : 
Gaspari et al. 20111 ^ 



star formation (jSchawinski et al.l 120071: iRoseboom et al.l 
I2009L but see also iShin et al.l I2011D . Far-IR ob- 
servations have also found powerful, p lausibly AGN 
drive n outflows in 0H~ and CO 1-0 (F eruglio et al.l 
I20TQI : iChung et al.l I20TTI : [S turm et al. 20li|) in local 
ULIRGs and QSOs that could exhaust the fuel sup- 
ply for star formation within ^10^ years, while X- 
ray observations have found mildly relativistic outflows 
in radio-quiet and radio- loud AGN, with a n origin 
close to the SMBH fChartas e t al.l ^03; Brait o et all 
120071 : fTombesi et al.llSoiOf ). Further observations have 
provided indirect e vidence for ^quasa r mode' feedback 
galaxies (iDunn et al.l 120101: lAlexander et al 

mi 



m some 



20101 iMuller-Sanchez et al.l 
2OIII ) , a nd for rad io -mo de 
clusters (iBest et al. 2006; 



2OIII : iRupke fc Veilleux 



feedback in galaxie s and 
Mc Namara fc Nulse nI 120071: 



Observationally, evidence for AGN feedback is mount- 
ing. Studies of early type galaxies suggest that 
some form of feedback may have terminated their 



Mittal et al.l 12009; Nes vadba et al.l 12010: Werner et al 

1 I20TTI : iMa et al 



2010; Ehlert et al. 201lt iShabala et al 
i201L : .Werner et al.i i201l|). Controversies do 



remam. 



^^.^ , however, 

For example. iJahnke fc Macci6l ()2011f ) suggest 
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Fig. 1. — Observed- frame fits to the optical through MIPS SEDs for objects 1-6 in Table [T] The solid line is the combin ed b est-fit model, 
while the dashed and dotted lines are the starburst and AGN components, respectively. A description of the models is in ^3.11 The number 
in the top right hand corner of each plot is the ID number in Tabled 



that there is no need for A GN feedback to explai n the 
Bulge-SMBH mass relation, I Amnions et alJ (|2Qllf ) sug- 
gest that only a small fraction of intermediate luminosity 
AGN at z ^ 1 can be undergoing galaxy- wide AGN feed- 
back, and iLutz et all ()2QlQf ) find that there is no simple 
inverse relation between star formation rate and AGN 
luminosity in X-ray selected AGN. It is also worth men- 
tioning that the eflFectiveness of AGN feedback may de- 
pend on the mode of star formation, since spatially com- 
pact, 'bulge' star formation is likely easier to turn off 
than spatially extended, 'disk' star formation. 

Our group has been examining the role of AGN feed- 
back by looking at systems in which such feedback may 
be ongoing. To do so, we have been studying the 



'FeLoBAL' class of QSOs ([Hazard et al.lll987D . which 
comprise part of the Broad Absorption Line (BAL) QSO 
populatioiF^. We selected this population for three rea- 

BAL QSOs show broad, deep absorption troughs in their rest- 
frame UV spectra. They come in three subtypes. High Ioniza- 
tion BAL QSOs (HiBALs) show absorption in CIV A1549A, NV 
A1240A, SilV A1394A & Lya. Low Ionization BAL QSOs (LoB- 
ALs) additionally show absorption in Mgll A2799A and other lower 
ionization species. Finally, the FeLoBAL QSOs, in addition to 
showing all the absorption lines seen in LoBALs, also show weak 
absorption from any excited term of Fell (e.g. Fell* A2400,2600A), 
& also Fell A2750A, Felll A1895,1914,1926A. See e.g. ILvndsl 
1967; Wevmann et al. 1991; Green & Mathur 1996; BeckereOD 
1997 ; Schmidt Hines 1999 ; Hall et al. 2002 ; Lacv et al. 200^ : 
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sons. First, the UV absorption troughs are unambiguous 
signatures of radiatively driven outflows powered by an 
AGN. They have velocities of up to 0.2c, widths of at 
least a few thousand km s~^), and are usually very deep 
or black, implying high column densities moving at high 
velocities, w hich plausibly iniplies hi^h mass- loss rates 
(lArav et al.| |1994VProga et al.ll2QQQl: IChartas et aIll2QQ3l : 
iCrenshaw et al. 2003). Therefore, they cannot be driven 
by even the most extreme starbursts. Furthermore, in 
at least some FeLoBAL QSOs the out flows m ay extend 
up to several Kp c into the host galaxy (|Arav et al. 2008; 
iMoe et al.ll2009l : iBautista et al .1120101 : iDunn et al...2010i ). 

Trump et aP ^200^; 'Gallagher et alj [20071 : ICasebeer eFall [20081 : 
Gibson et al. 2009 : Leighlv et al. 2 01ll . 



Second, they ar e invariably redden ed objects with high 
IR luminosities (Farrah et al.""2007l, and sometimes har- 
bour intense starbursts (Farrah et al. 2010), suggesting 
they may be a transition phase in the lifetime of an AGN 
(FeLoB AL features ar e also occasionally seen in ULIRGs, 
see e.g. [Farrah et al.| [2005). Third, recent results have 
shown that FeLoBAL QSOs are much more common at 
z ^ 0.5 t han was originally t hought, by up to a factor 
often (|Dai et al.[ [20081 : lUrrutia et al.[ [20091 : [Allen et al.[ 

[loiiD. 

In this paper, we compare the strengths of the 
AGN-driven outflows to the luminosities of obscured 
star formation in a sample of 31 FeLoBAL QSOs, se- 
lected purely on the basis of their rest-frame UV spec- 
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Fig. 3. — Observed-frame fits to the optical through MIPS SEDs for objects 13-18 in Tabled Details are the same as for figure ^ 



tral properties. We combine da ta from the Spitzer 
space telescope (Wer ner et alJ 120041) with d ata from the 
Sloan Digital Sky Survey (SDSS, lYork et al. 20001 the 
Two Micron All-Sky Survey f2MASS. ISkrutskie et all 
2006h. the UKIR T Infrared Deep Sky Survey (UKIDSS, 
Lawrence et al.l ^007), and from the Wide-Fie l d In- 
frared Surv ey Ex plorer (WISE, "Wright et alj [201(1 
Ijarrett et al.ll2QTTI ) in order to measure the luminosities 
of both obscured AGN activity and star formation. We 
measure the strengths of the AGN-driven outflows from 
the UV absorption troughs in the SDSS spectra. We as- 
sume a spatially flat cosmology, with Hq = 70 km s~^ 
Mpc"\ ft = 1, and = 0.7. We use the term "IR 
luminosity" to refer to the luminosity integrated over 1- 



1000/im in the rest-frame. We quote luminosities in units 
of bolometric Solar luminosities, where = 3.826 x 10^^ 
Watts. 

2. ANALYSIS 

2.1. Sample Selection 

We aimed to select a sample of FeLoBAL QSOs 
purely on the basis of their rest-frame UV spectral fea- 
tures. Furthermore, since measuring BAL properties 
is not straig:htforward (see 34.21 and iHall et al.l [20021 : 
[Trump et alT[2006l : [Gibson et al.[[2009[ : [Allen et al.[[20Tl[ ) 
we required the sample to have optical spectra from the 
same source, and BAL measurements already in the lit- 
erature. 
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Fig. 4. — Observed-frame fits to the optical through MIPS SEDs for objects 19-24 in Table [T] Details are the same as for figure [T] 



Accordingly, we chose our sample fro ni the SPSS. 
We started with the six SDSS objects in iFarrah et al.l 
^007 (hereafter F07), excluding ISO J005645. 1-273816 
and LBQS 0059-2735 as neither lie within the SDSS 
survey. These six objects were originally selected with 
the only constraint being that their redshifts satisfied 
1.0 < z < 1 8. We then s electe d a further 25 FeLoBAL 
QSOs from iTrump et aP (|2006l ). We imposed the same 
upper redshift cut of z = 1.8 to ensure that observed- 
frame 160/im remains close to the peak of the far-IR 
emission for the (monolithic) dust temperatures of > 50K 
expected within ~ lOOpc of an AGN. We slightly reduced 
the lower redshift cut from the F07 sample to z = 0.8 to 
ensure that the SDSS spectra always contain the Mg II 



BAL (see ^4.2.ip . Hence, the final sample comprises 31 
objects. This sample should be a random subset of the 
FeLoBAL QSO population between 0.8 < z < 1.8, with 
no selection on IR luminosity. The sample is listed in 
Table HI 

We are however not certain that our sample are all 
FeLoBAL QSOs, sinc e BAL features are difiicu lt to iden- 
tify (see go and e.g. lAppenzeller et al.l[2005l ). Two ob- 
jects in particular are problematic; SDSS J033810 and 
SDSS J233646. SDSS J033810 is discussed in detail by 
iHah et al.ll200l in summary, it is not certain whether 
the BALs in this object are real, or a phantom of a pe- 
culiarly shaped continuum. SDSS J233646 on the other 
hand is clearly an FeLoBAL QSO, but is in a binary sys- 
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Fig. 5. — Observed- frame fits to the optical through MIPS SEDs for objects 25-30 in Table [T] Details are the same as for figure [T] 



tern with a separation of < 2^' (jGregg et al.ll2QQ2l ). We 
have assumed that the IR emission comes solely from 
the FeLoBAL QSO (object B in Gregg et al), but it 
is probable, particularly at 24/im, that there is con- 
tamination from the companion. Two further objects 
are worth mentioning; SDSS J081312 & SDSS J105748. 
These two objects show absorption from Fell but not 
Fell*; the feature originally identified as FeII*A2600A 
was subsequently shown to be a blend of FeIIA2580A 
MnIIA2577,2595,2606A and narrow MgIIA2799A with 
an uncertain relation to the QSO absorption. Since the 
original definition of an FeLoBAL Q SO requires the de- 
tection of Fell* (Haza rd et al.l [19871 ), these two objects 
would not be formally classified as FeLoBAL QSOs if ex- 



treme strictness was imposed. They do however both 
show Fell absorption, so we treat them as FeLoBAL 
QSOs. 

2.2. Observations 

A description of the observations of the six objects from 
F07 can be found in that paper. The 25 new objects were 
observed in cycle 5 with t he Multiband Im aging Pho- 
tometer for Spitzer fMIPS. [Rieke et all [2004 ) at 24/im, 
70/im and 160/im. The small field size was used for all 
three channels, using the default pixel scale at 70/im. At 
70/im and 160/im we chose the same exposure parameters 
as used in F07; seven 10 second cycles at 70/im and four 
10 second cycles at 160/im. At 24/im we used a shallower 
total exposure time of two 10 second cycles per source. 
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Fig. 6. — Observed-frame fit to the optical through MIPS SEDs 
for object 31 in Table [T] Details are the same as for figure [T] 

sufficient to detect, at > 5cr, a 0.5mJy source, since all 
the objects in F07 were strongly detected at 24/im. 

We reduced the data using the MOPEX software pro- 
vided by the Spitzer Science Centre, which performs 
standard tasks such as image co-addition, sky and dark 
subtraction, and bias removal. We used the default 
MOPEX parameters for MIPS small field observations. 
To check the quality of our reduction we compared it to 
that from the automated pipeline and found them to be 
consistent to within a few percent. We measured fluxes 
in all three channels using two methods. For the obser- 
vations in which the source was clearly detected - all the 
24/im data and about half of the 70/im data - we used the 
APEX package within MOPEX with the default param- 
eters and PRFs for point-source photometry. For the ob- 
servations in which the source was weakly or not detected 
we found that PRE fitting photometry was not suitable, 
so w e used the Sextractor package (Bertin & Arnoutsl 
Il996f ) in aperture photometry mode, with point-source 
aperture corrections provided by the Spitzer Science Cen- 
ter. For undetected sources we measured the flux within 
the aperture, and then subtracted the average back- 
ground flux. We compared the results from APEX PRF 
fitting photometry and Sextractor aperture photometry 
for the brighter sources, and found that they agreed to 
within a few percent, with no obvious systematic differ- 
ences. We therefore use the resulting mix of aperture 
and PRF fluxes in our analysis. 

3. RESULTS 

The MIPS photometry is presented in Table O We 
combine the MIPS data with archival photometry from 
SDSS, 2MASS (or UKIDSS if available), and WISE. This 
archival data is presented in Table [TJ Additi onal sub- mm 
photo metry for two objects is presented in iLewis et al.l 

As WISE returned data only recently, we provide more 
detail on the WISE data. WISE completed its first full 
coverage of the sky in July 2010. A Preliminary Re- 
lease Catalogs covering 57% of the sky, was made avail- 
able to the community starting April 14, 2011. Point 
sources were extracted with 5-cr sensitivities greater than 
0.08, 0.11, 0.8, and 4 mJy, respectively for the four 
bands. Photometric calibration, source counts, colors 

http:/ /wise2. ipac.caltech.edu/docs/release/prelim/ 



and popul ation statistic s at t he ecliptic poles are pre- 
sented bv iJarrett et al.l ()2011f ). Of our 31 sources, ten 
come from the Preliminary Release Catalog, and the 
rest come from the proprietary WISE First Pass Internal 
Source Database, which performs internal verification, 
quality and photometric analysis. In the cases where 
there were multiple measurements of the same source 
(due to the overlap between WISE orbit-to-orbit scans) 
we used quality indicators to choose the best measure- 
ment. Overall, 25 of our sources have > 2a detections in 
at least one of the four WISE bands. 

We measure total IR luminosities, and the contribu- 
tions to the total IR luminosities from star formation 
and AGN activity, by fitting the optical through MIPS 
photometry for each object with radiative transfer mod- 
els for AGN and starbursts, following the methods in 
Farrah et al."2003' and F07. We describe the models and 
the fitting methods in g3l]& 

3.1. The Models 

For the AGN ni o dels, we follow 
lEfstathiou fc Rowan- RobinsonI ()199 5') and use a ta- 
pered disk dust distribution, in which the disk thickness 
increases linearly with distance from the central source 
in the inner part of the disk but tapers off to a constant 
height in the outer part. The tapered disc models 
include a distribution of grain species and sizes, multiple 
scattering and a density distribution that follows 
where r is the distance from the central source. The 
models assumed a smooth distribution of dust, so they 
are a good approximation of the density distribution in 
the torus if the mean distance between clouds is small 
compared with the size of the torus. These models have 
been successful in fitting the SEDs of several classe s of 
AGN, and ULIRG s ([Alexander et al. 1999; Ruiz etld] 
12001; Farrah et al. 2002 ;IVerma et al.ll2002UFarrah et al.l 
12003: Efstathiou fc Siebenmorgenll2005f ). In this paper 
we use a grid of models with four discrete values for 
the equatorial lOOOA optical depth (500, 750, 1000, 
1250), three values for the ratio of outer to inner disc 
radii (20, 60, 100) and three values for the opening 
angle of the disc (30^, 45^ and 60^). The spectra are 
computed for 74 inclinations which are equally spaced 
in the range to 7r/2. For comparison, other work 
on radiative transfer modeling of the d ust dist ri butio n 
in AGN has been pr esente d by Pier fc KrolikI (11992[) : 

[Granato fc Danesd (ll^gl ): iNen kova et al.l ( 2002li : 

Dullemond fc van Bemmell (j2005f ): iHonig et all (j2006l i: 
>Schartmann et al.l (j2008f ). 

For the starb urst models, we combine the 

iBruzual fc Charlotl (|2003f ) stellar population syn- 
thesis model with a radiative transfer code that includes 
the effect of s mall grains and PAHs. the u pdated dust 
model of Efsta tliiou fc SiebenmorgenI (|2009l ). and a sim- 
ple evolutionary scheme for the molecular clouds that 
constitute the starburst. The model predicts the spectral 
energy distributions of starburst galaxies from the ul- 
traviolet to the millimetre as a function of the age of the 
starburst and the initial optical depth of the molecular 
clouds. These models have been successfully used to 
fit th e SEDs of several classes of star-forming galaxy 
(e. g. [Farrah et aLl 120031 ), as well as the Tor k^ diagram 
of iSpoon et al.l l2007l ([Ro wan- Ro binsonfc Efstathioul 
[lOOi). For comparison, other starburst models have 
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Fig. 7. — An example of the method used to derive the IR luminosities for each object. Here we show the data for one object, SDSS 
J221511. 93-004549. 9. The left panel shows the Xred values and resulting luminosities for all possible combinations of starburst and AGN 
template fits. The green points are the total IR luminosities, while the red and blue points are the AGN and starburst luminosities, 
respectively. The right panel shows the differential probability distribution functions (for total, AGN and starburst luminosity) that arise 
from the Xred weighted combination of all the solutions in the left panel. 



been developed by Rowan-Robinson & Crawford" 
( 1989[) : iKruegel fc Sieben morgen (1994); Silva et al. 
(119981): rTakagj et all l2QQ3f) : IDopita et all (|2005|); 
iSiebenmorgen fc Kriigell (|2QQ7l ). 

3.2. Fitting & Results 

We fit all possible combinations of starburst plus AGN 
models to the SDSS through MIPS photometry of each 
object in order to extract starburst and AGN luminosi- 
ties. 

The models do not include a prescription for the broad 
absorption features in the rest-frame UV that arise from 
radiatively driven outflows, so we do not include the 
SDSS spectra in the fits. We do not include photome- 
try shortward of '^0.35/im in the rest-frame for the same 
reason; this means that we do not use the Ugri pho- 
tometry for any object, and in some cases do not use 
the z band data either. We do however use photome- 
try at rest-frame wavelengths of 0.35 < A(/im)< 1; while 
we are interested in starburst and AGN luminosities at 
rest-frame 1-1000/im, including these data limits which 
models are acceptable fits in the optical, and thus limits 
which models can be used to fit the IR data. For the 
MIPS data we fit to the measurements in all cases even 
if they are of low significance, or negative. We allowed 
the contribution from each component to vary from 0% 
to 100% to see if the IR emission was consistent with 
arising from a single origin. The model libraries span 
several free parameters ( §3.ip . but our data are not com- 
prehensive enough to constrain all of them. Therefore, 
we use the complete model libraries to determine only 
the range in total, starburst and AGN luminosities that 
are consistent with the data. 

The best-fit SEDs are presented in Figures 1-6. In most 
cases the models provide an excellent fit to the data. We 
obtain xled ^ ^ 12/31 objects, and 1 < xled ^ ^ 
in 15/31 objects. The remaining four objects have 2 < 
Xred ^ ^- most difficult points to fit are usually the 
near-IR and/or the two shorter wavelength WISE bands, 
which are somewhat underpredicted in many cases. This 
could be due to the still preliminary WISE photometric 



calibration, and/or the lack of a host galaxy contribution 
in our models, so we defer exploration of this until more 
comprehensive near-IR data are available, and the WISE 
calibration is refined. In nearly all cases the > 24/im 
photometry is well fitted by the models, though in a few 
objects the fit is relatively poor at these wavelengths. 
This effect is most noticeable in objects 3, 11, and 24, in 
which the 70/im flux is significantly underpredicted. This 
could be due to an additional 'hot' dust component or an 
extremely strong spectral feature that is not reproduced 
in the models. Since the models still reproduce the slope 
and approximate normalization of the far-IR SEDs even 
for these objects, we consider the results to still be usable. 

To extract the most probable total, AGN and starburst 
luminosities for each object, and their confidence inter- 
vals, we first construct a discrete probability distribution 
function (PDF) for the total, AGN and starburst lumi- 
nosities of each source. We weight each fits contribution 
to the PDF by its xled value. From these PDFs, we 
then construct cumulative distribution functions, from 
which the most probable luminosities, and their confi- 
dence intervals, can be derived straightforwardly. Ide- 
ally, we would show the PDFs for all objects. This would 
however be unwieldy, so, as an example of the method, 
we show the distribution of Xred values and the total, 
starburst and AGN PDFs for one object in Figure [71 
In most cases the PDFs have a single peak and an ap- 
proximately gaussian shape. In a few cases though there 
are minor secondary maxima, and/or the shape is non- 
gaussian. We therefore quote as the positive and negative 
errors the 90% confidence intervals. The luminosities, 
their confidence ranges and the Xred values of the best 
fits are presented in Tabled 

4. DISCUSSION 

4.1. Infrared Luminosities 

We find that FeLoBAL QSOs are luminous in the IR. 
All of our sample have (best-fit) total IR luminosities 
(LTot) in excess of lO^^L©, with nine objects, or 29%, 
exceeding lO^^L©. There is no other sample to which 
we can directly compare ours, since no sample matched 
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Fig. 8. — Left: Total IR luminosity vs. contribution to the total IR luminosity from the starburst. Right: AGN luminosity vs. starburst 
contribution to the total IR luminosity. In neither plot do we see a clear correlation, though it is interesting that the starburst contribution 
drops at L^gn ^ lO^^L©. 



to ours in optical continuum luminosity and redshift has 
been observed in the mid/far-IR. We do however find 
that FeLoBAL QSOs are mor e IR- lumino us than the 
Palomar-Green (PC) QSOs in i Haas et aP (|2QQ3l In- 
stead, they appear to have comparable IR luminosities 
to th e wider population of BAL QSOs (jGa llagher et alJ 
I2QQ7D . reddened QSOs (jGeorgakakis et a L 2009), the 
most IR-luminous members of the general QSO pop- 
ulation (Lutz et al. 2008; Orellana et al. 2011), and to 
ULIR Gs (Genzel et all 119981 : IFarrah et all 120021 [20031 : 
iDesai et aL 2007). 

The dominant power source behind the IR emission 
is, in most cases, AGN activity. A pure AGN is either 
the most likely power source, or consistent within the 
90% confidence interval, for 35% of the sample (11/31 
objects). A starburst component is required (at > 90% 
confidence) for the remaining objects, but in only twelve 
of these objects is the starburst more luminous than 
lO^^L©, and in only three objects is the starburst more 
luminous than the AGN. The mean AGN contribution 
to the total IR luminosity is ~ 76%, comparable to that 
seen in local ULIRGs with ^wa rm' IR colour but lower 
than that seen in PG QSOs (jVeilleux et al.ll2009f ). The 
spread in AGN contribution to Lrot is however wide, 
spanning 0.2 to 1.0. We find a strong correlation be- 
tween LTot and the AGN IR luminosity (Lagn)^ with a 
Spearman rank correlation cofficienlF^ of p = 0.92 and 
a significance of deviation from zero of P < 0.001. We 
also find a correlation between L^ot and starburst IR 
luminosity (L^^), with p = 0.66 and P < 0.001. Con- 
versely, we find no trend between L^ot and the starburst 
contribution to L^ot ifss)^ with p = —0.23, P = 0.21 
(Figure [8] left). Finally, we find a weak anticorrelation 
between Lagn and fsB {p = -0.47, P = 0.007, Figure 
[8] right). This is most apparent when the AGN IR lumi- 
nosity exceeds ~ lO^^L©, at which point the starburst 

That is, those objects where the IRAS 12/im to IRAS 60/im 
flux ratio is greater than 0.2 - this is usually interpreted as an AGN 
supplying at least a plurality of the IR emission. 

We also evaluated all the correlations presented here using 
the Kendall tau rank correlation coeflicient, and obtained similar 
results in all cases. 



contribution decreases noticeably. 

The star formation rates for the objects with starbursts 
detected at > 90% confidence, computed by extract- 
ing the monochromatic 60pm lu minosities from the SEP 
fits an d then using equation 7 of iRo wan- Robinson et al.l 
(|1997[ ). lie between several hundred to a few thousand 
solar masses per year. F or the six objects in common 
with "Farrah et al.l (|2010l ). these star formation rate es- 
timates agree with those from the IRS spectra, though 
the large systematics on both measurements renders this 
comparison of little value. 

We now examine the IR emission from AGN and star- 
bursts in FeLoBAL QSOs using a different approach. 
We use all the solutions from all the individual fits to 
construct a single PDF for the whole samplJ^ (Figure 
[9]). From this, we reach similar conclusions to those de- 
scribed above. An AGN supplies the bulk of the IR emis- 
sion in the majority of cases. Starburst activity on the 
other hand is significantly less luminous; we find that 
the starburst is fainter than lO^^L© in just over half 
{^6.7^.1'^%) of cases, and almost never exceeds lO^^L© 
(Table il first row). 

Overall, we find compelling evidence from the IR prop- 
erties of FeLoBAL QSOs that they are, as a class, dis- 
tinct from the general QSO population. FeLoBAL QSOs 
have higher IR luminosities and a smaller average AGN 
fraction. Instead, the IR properties of FeLoBAL QSOs 
are consistent with those of the reddened & LoBAL QSO 
populations. 

The relationship between FeLoBAL QSOs and 
ULIRGs is harder to discern. It has been suggested 
(F07) that FeLoBAL QSOs are, as a class, a ULIRG- 
to-QSO transition phase. There is one result from our 
study that supports this, namely that the mean AGN 
fractional IR luminosity of FeLoBAL QSOs is similar to 
that of 'warm' ULIRGs. Conversely, if FeLoBAL QSOs 
were such a transition phase then we may have expected 
a higher fraction of them to host luminous starbursts, 
up to ~50% if they represented the entirety of the star- 

This is more appropriate than taking the individual PDFs and 
multiplying them together to make a single PDF, since the quality 
of the best flts from object to object varies signiflcantly 
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Fig. 9. — Differential PDF for the whole sample, constructed from all fits to all objects, weighted by their individual Xjied ^^ilues. The 
PDF for the total IR luminosity is plotted in green, for the AGN luminosity in red, and for the starburst luminosity in blue. The top panel 
shows all three PDFs overlaid, while the bottom row shows the individual PDFs and their errors. The errors were estimated via jack-knife 
resampling, removing one source at a time and deriving the resulting error estimates, but were performed for the bins individually. Since 
the bins are correlated, these errors are an overestimate of the overall uncertainties in the PDFs. 



burst to AGN transition (see e.g. iFarrah et al.ll2QQ9l ). 
Now, this fraction is an upper hmit to what we might 
see, and such a high fraction with luminous starbursts 
is (just) consistent with the PDF in Figure [9l Other 
evidence for this comes from the unusual mid-IR spec- 
tral shapes of FeLo BAL QSOs i n com parison to other 
IR-luminous QSOs IFarrah et aP (\201(f ). From the re- 
sults in this paper however, the scenarios that FeLoBAL 
QSOs are (a) randomly drawn from the reddened QSO 
population, and (b) a ULIRG to QSO transition phase, 
are both plausible. Since the former scenario is simpler, 
we conclude that, while some fraction of reddened/IR- 



luminous QSOs are almost certainly the endpoint of a 
ULIRG-QSO transition, we see no evidence in this paper 
that demands that FeLoBAL QSOs are more likely to 
be such a transition stage than the (presumably) parent 
reddened QSO population. 

4.2. Starbursts & Outflows 

We now examine whether or not there is a relationship 
between the AGN-driven outflows and the obscured star 
formation in our sample. In §4.2.11 we discuss how to 
measure the strength of the BAL features. In §4.2.21 we 
examine the relationship between AGN-driven outflows 
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Fig. 10. — Absorption strength vs (top left) total IR luminosity, (top right) AGN luminosity, (bottom left) starburst luminosity, and 
(bottom right) starburst contribution to the total IR luminosity. Objects with a starburst detected at > 90% confidence are plotted in 
colour; in blue if they exceed IO^^Lq, and in green if they are fainter than this luminosity. Objects where the 90% confidence interval on 
the starburst luminosity encompasses zero are plotted in black. 



and obscured star formation using the best-fit luminosi- 
ties and their confidence ranges to start with, and then 
examine it using all the information in the PDFs. We 
interpret our results in §4.2.31 

There is one point we note first. We cannot here mea- 
sure how much kinetic energy the outflow is injecting into 
the ISM of its host galaxy. Even with well resolved BALs 
in multiple species, such a measurement involv es in- 
depth radiative transfer calculations (e.g. Casebeer et al.l 
I2OO8). The SDSS optical spectra used here are how- 
ever of low resolution, and usually contain BALs only 
for one or two species. So, even attempting such calcu- 
lations for our sample is futile. It is however plausible 
that the depth, width, and velocity offset of BALs in 
a single species, as quantified in mea sures such as Bal- 
nicity Index and Absorption Index (H azard et al.l 119871 : 
iHah et al.ll2QQ2l : lTrump et al.fe OOG) do scale with increas- 
ing outflow strength. So, we here use the properties of 
the Mg II BAL features solely as an estimate of relative 
outflow strength within the sample. 

4.2.1. Measuring Outflow Strengths 

The strength of absorption of the BALs in a given 
species is traditionally defined as the total velocity width 



over which the absorption exceeds a minimum value. 
This however gives rise to four problems with interpreta- 
tion. First, different authors use different parametriza- 
tions for this strength, which in some cases can determine 
whether an object is identified as a BAL QSO (of any va- 
riety) or not. Second, the derived depths and widths of 
the BAL troughs are sensitive to the choice of contin- 
uum level, which can vary significantly between auto- 
mated and manual measurements ev en for the same in- 
put data (e.g . conipare the results inlTrump et ani2QQ6l : 
iGibson et al][2QQl lAlIen et al.ll2QTT h. Third, since our 
sample is at 2; < 1.8 we see absorption only in a limited 
number of species. Fourth, the BALs in some QSOs show 
temporal variatio n on approximately decadal timescales 
(jArav et al.ll2QQlf ). 

The latter two issues cannot be addressed with the 
data we have. To mitigate the first two as far as is pos- 
sible we proceeded as follows. First, we used the follow- 
ing parametrization for the Absorption Strength (AS) of 
BALs: 



AS-- 





\l ^^""^ 


Bdv 


J Vo 
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TABLE 3 

Examples of different absorption strength 
measures from equation [u ordered by decreasing 
strictness. 



Name 


vo 


VI 


a 


Min. Velocity Width 


ASo (BI)- 


3000 


25000 


0.9 


2000 


ASi 


3000 


29000 


0.9 


2000 


AS2 (BIo)^ 





25000 


0.9 


2000 


AS3 





29000 


0.9 


1000 


AS4 (AI)^ 





29000 


1.0 


1000 



Note. — The names in parentheses are those of measures 
used in the Hterature. We do not use ASi or AS3 here, but 
present them to illustrate other possible parametrizations. 
^ Balnicity Index (IWevmann et al.lfTOQll : IGibson et al.|[2QQ9l : 
I Allen ena1f2QTlh . 

^ Modified Balnicity Index (iGTbson et al.|[20Q9l). 

^ Absorption Index JHall et al.ll20Q2l : iTrump et al.ll2QQ6l) . 

where v is velocity, f{v) is the (normahzed) flux at that 
velocity, and a is a scaling factor. The quantity B is set to 
unity if the absorption is below 10% of the continuum and 
if the width of the trough is greater than some minimum 
value, otherwise it is set to zero. We deflne the systemic 
redshift to be at zero velocity, and positive velocity to be 
blueward of this redshift. Equation [T] includes all previ- 
ous deflnitions of BAL strength as special cases (Table [3j). 
We chose the AS2 parametrization as our primary mea- 
sure since it is a reasonable compromise between strict- 
ness and inclusiveness. Second, we measure absorption 
strengths for the same, single species across the whole 
sample. We chose MgIIA2799A as it is the only species 
present in all the SDSS spectra of our sample, though 
we note that it can be contami nated by Fe II absorption 
(see Table 1 of iHall et al1l2QQ2i) . We remeasured ah the 
Mgll absorp tion strengths by ha nd, using the methods 
described in lUrrutia et al.l (|2QQ9[ ). Our measurements, 
together^with one se t of comparison measurements from 
iGibson et al.l (j2QQ9f ). are given in Table El We check 
our absorption measures against those in the literature, 
and explore the effects on our results of using different 
measurements of absorption strength, in the appendix. 
The formal error on the absorption strengths from the 
flts is usually of order 100 km s~^, but there are sig- 
niflcantly larger systematic unertainties (see appendix) 
that are hard to quantify. These systematic uncertainties 
should however not dramatically change the absorption 
strengths of the sample relative to each other, as long 
as the measurements are done in an internally consistent 
way. 

4.2.2. Outflows vs. Infrared Properties 

We flrst use the luminosities in Table [2] to investigate 
if absorption strength depends on total, AGN and star- 
burst IR luminosity. We flnd no correlation between ab- 
sorption strength and LTot {p = 0.31, P = 0.09, Fig- 
ure [To] top left), or between absorption strength and 
^SB {p = -0.10, P = 0.58, Figure [JO] bottom left), 
though we do flnd a hint of a correlation between ab- 
sorption strength and Lagn {p = 0.39, P = 0.04, Fig- 
ure [10] top right). It is also interesting that (a) the dis- 
persion in absorption strength is greater at higher total 
and AGN luminosities, and (b) all but two of the lumi- 
nous (> lO^^L©) starbursts lie at AS2 < 5000 km s~^. 
We conclude that the mechanism that determines the 



strength of the outflows is not directly responsible for 
heating the dust near the AGN, and does not have a 
strong effect on the absolute luminosity of the starburst. 

If however we plot absorption strength against fsB 
then we see a weak but clear anticorrelation (p = —0.49, 
P = 0.005, Figure [TQ] bottom right). Moreover, there 
seems to be a change in the distribution of starburst 
contributions with absorption strength at AS2 — 3500km 
s~^; all the systems with AS2 > 3500km s~^ have a star- 
burst contribution of less than 25%, while the systems 
with AS2 < 3500km s~^ have a wide dispersion in star- 
burst contributions, from 0% to ~ 80%. 

To estimate the probability that the distribution of 
starburst contributions changes at AS2 = 3500km s~^ 
we employ the two-sided Kolmogorov-Smirnov test. We 
flnd the distributions of the objects above and below 
AS2 = 3500km s~^ in the bottom right panel of Fig- 
ure [10] are different at 99.84% confldenc43. The number 
of objects in the two subsamples is however low enough 
that the Kolmogorov-Smirnov test can be unreliable. So, 
we employ a cruder test. We take the null hypothesis 
to be that there is no correlation between Absorption 
Strength and fsB^ and that the systems with fsB < 25% 
represent the 'true' underlying distribution. The prob- 
ability of flnding all eight systems with fsB > 25% at 
AS2 < 3500 km s-i is then (8/19)^) 1%, i.e. a ::^99.0% 
probability of a difference. We regard this latter flgure 
as a more reliable measure of the signiflcance of a differ- 
ence. Hence our flnding all of the fsB > 25% systems 
at AS2 < 3500 km s~^ is only weak evidence that an 
anticorrelation exists. 

Nevertheless, this result is consistent with the idea that 
the radiatively driven outflows negatively affect star for- 
mation. The systems with fsB < 25% would then be 
those in which an outflow has curtailed star formation, 
and those in which such an outflow has subsequently 
waned, making the observed dispersion in absorption 
strengths wide. The systems with fsB > 25% would 
be those in which an outflow has started to develop, but 
has not yet affected the starburst. 

To explore the relationship between absorption 
strength and infrared properties further, we use all of 
the information in the PDFs, in a manner similar to that 
used in §4.11 We adopt a boundary value of absorption 
strength motivated by Figure [TO] of AS2 = 3500km s~^, 
and divide our sample into two subsamples at this bound- 
ary. 

We flrst construct PDFs of Lagn and L^^ for these 
two subsamples, and extract from them the probabilities 
of obtaining luminosities in excess of certain values (rows 
2 & 3 of Table [4|) . We see similar results to those seen in 
the flrst three panels of Figure [10] We see no convincing 
differences in the starburst luminosities between the low 
and high absorption strength subsamples, compared to 
either each other or to the sample as a whole. We also 
see that we are only marginally more likely to see AGN 
with Ljji < in the low absorption strength sub- 

sample. 

If however we consider the PDFs for the contribution 
of the starburst to the total IR luminosity for the two 
subsamples (Figure [TT]) . then we see a clear difference. 

Employing the conceptually similar Mann- Whitney test gives 
a comparable result 
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TABLE 4 

Probabilities of obtaining Starburst & AGN luminosities below certain 

BOUNDARIES. 



Selection 


< IO" 


Probability of Obtaining 

< 1012 < 1012.5 < 1012 


(Lo) 

< 1012-5 


All Objects 


28.3^^;^% 


56.7i^-j% 


83.7i^;^% 


21.4l?o%% 


57.0l^-J% 


AS2 < 3500 km s-i 
AS2 > 3500 km s-i 


27.7^^-^% 
29.6±^;2% 


55.21^-^% 
59.7±iio°% 


80.0+^-^% 


32.1+^4^1% 
0%' 


68.7+|-^% 
33.8±}0;^% 


Lagat < 1O12-5L0 
Lagn > 1O12-5L0 


38.8t^;^% 
16.2±2-J% 


69.4t^-^% 
43.21^0^^% 


93.7t^;g% 
69.8±^;2% 


37.3t{5^,% 
0% 


90.9+^-^% 
10.4±2;o% 



Note. — We give probabilities for the complete sample, the sample divided by absorption 
strength, and the sample divided by AGN luminosity (see also Table [61 . Errors were derived 
using jack-knife resampling, removing one source at a time and computing the ~ la confidence 
interval from all the resulting realizations. The subsamples divided by Lagn were divided on 
their peak luminosities, hence the non-zero probabilities of obtaining luminosities outside the 
boundaries. 
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Fig. 11. — Probability Distribution Function for starburst contribution divided according to their absorption strengths. The black line 
is the PDF for the whole sample, in red for the objects with AS2 < 3500 km s~i, and in orange for the objects with AS2 > 3500 km s~i. 
The orange line shows a much higher chance of obtaining a starburst fractional luminosity in excess of 25% than the red or black lines. 
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TABLE 5 

Probabilities of obtaining Starburst 
contributions to the total ir luminosity, 
above two boundaries. 



Selection 


> 25% > 50% 


All Objects 


50.3l^-^% 




AS2 < 3500 km 
AS2 > 3500 km 


67.3^^-^% 


26.6tl[i% 
< 1.5% 


Lagn < 1O12-5L0 
Lagn > 1O12-5L0 


57.8^^-1% 
38.7l?o'o% 


i8.6ti:i% 

15.1+^-0% 



Note. — As with Table ^ we give probabilities 
for the complete sample, the sample divided by ab- 
sorption strength, and our sample divided by AGN 
luminosity (see also Table [6]). Errors were derived 
using jack-knife resampling. 



The low absorption strength subsample shows a higher 
chance of a higher fsB than the sample as a whole, and 
a much higher chance than the high absorption strength 
subsample. We quantify this by extracting the proba- 
bilities of obtaining starburst contributions to the total 
IR luminosity in excess of 25% and 50% from the whole 
sample, and the two absorption strength subsamples (Ta- 
ble [5l rows 1-3). We find, at > 5cr significance, a higher 
chance of obtaining fsB > 25% and fsB > 50% in the 
low absorption strength sample compared to the high ab- 
sorption strength sample. These results do not change if 
we exclude the objects noted as potential contaminants 
in §2.H or if we exclude the objects with xjied > 2 in 
Table El 

4.2.3. Interpretation 

The anticorrelation that we observe between absorp- 
tion strength and contribution from star formation to 
the total IR emission is straightforwardly interpreted as 
the outflow from the AGN curtailing star formation in 
the host galaxy. There are however four other ways that 
we could see this anticorrelation. 

The first is that stronger outflows reflect an increase 
in the IR emission from the AGN, but have no effect on 
the starburst; if this is the case we would see a decline in 
fs B 5 but without there being a direct relationship behind 
the two phenomena. This possibility is apposite if there 
is a conspiracy of fortuitous timescales, where the peak 
starburst luminosity precedes the peak AGN luminosity 
by a few Myr. The second is that starburst activity sup- 
presses AGN outflows, SO when the starburst wanes (via 
a cause unrelated to the AGN) an AGN driven outflow 
can appear. The third is an observation bias; e.g. if a 
high fsB meant that the Mg II troughs were observed to 
be weaker than they really are. The fourth is a selection 
bias, e.g. if QSOs with strong BALs and strong star- 
bursts drop out of the in itial SPSS Q S O sel ection and 
were thus not included in iTrump et al.l ()2006f ) . 

We first consider the last three of these alternatives. 
The second alternative has traditionally been considered 
unlikely, the argument being that for it to happen then 
the starburst would have to act near the origin of the out- 
flow, but AGN broad-line regions are around five orders 
of magnitude smaller than starburst regions. Neverthe- 
less, it is not completely implausible. One possible sce- 
nario is as follows. If the ISM was initially dense and the 



SMBH was initially small, then the SMBH may not be at 
first capable of powering outflows that extend ~Kpc into 
the host galaxy, but as the density of the ISM was re- 
duced by the ongoing starburst (which thus waned due to 
the reduction in fuel supply) and the SMBH grows, then 
large-scale outflows would subsequently appear. This is 
not a scenario we can test, but it would likely require a 
serendipitous conjunction of ISM and SMBH parameters, 
SO we do not consider it further. 

The third possibility is one that we again cannot test, 
so we cannot formally discount the possibility of a very 
large population of OB stars acting to suppress the ob- 
served depth of the UV absorption troughs. Conversely, 
the rest-frame UV spectra of starbursts in ULIRGs re- 
veal continua that are at least an order of magnitude too 
weak to provide this effe ct, and sometimes show absorp- 
tion in the same species (Farrah et al.""2005"). The fourth 
possibility is also not testable, but the SDSS is now turn- 
ing up FeLoBAL QSOs in large numbers, and the initial 
QSO followup colour selections are fairly relaxed, so we 
do not consider this possibility likely either. 

The first alternative is however one that we can test, 
from which we propose that it is unlikely as well. The 
test is as follows. If it is the case that outflow strength 
is a proxy for AGN luminosity, then we should see a big- 
ger difference between the starburst contribution PDFs 
for subsamples divided by AGN luminosity than between 
subsamples divided by absorption strengtlF^. In Fig- 
ure [12] we show the starburst contribution PDFs for 
two subsamples divided by AGN luminosity, one for ob- 
jects with Lagn > and one for objects with 
Lagat < Qualitatively, the difference between 
the PDFs in Figure [12] is weaker than the difference be- 
tween the PDFs divided by absorption strength in Figure 
[TTl If we extract probabilities of obtaining the same star- 
burst & AGN luminosities, and starburst contributions 
as we did for the absorption strength subsets (4th & 5th 
rows of Tables |3]&[5|) then we see three interesting results. 
First, lower luminosity star formation is now more likely 
(at > 3<j) to be seen in the lower luminosity AGN sub- 
sample (e.g. forP(L55 < lO^^Lo): theLAGN < lO^^-^L© 
subsample is 69.4^4;^% while the Lagn > lO^^-^L© sub- 
sample is 43.2l^Q^4%). Second, we are more likely (albeit 
only at just over 2a) to obtain a smaller starburst con- 
tribution by selecting high absorption strength systems 
than we are by selecting high AGN luminosity systems 
(e.g. for F{fsB > 25%): the AS2 > 3500km s'^ subsam- 
ple is 17.8tg 5% while the Lagn > lO^^-^L© subsample 
is 38.7l^Q^Q%). In other words, we are more successful 
in finding systems with a large starburst contribution 
to the total IR emission by selecting on weak outflows 
than we are by selecting on low AGN luminosity. Third, 
the probabilities of a starburst contribution in excess of 
25% (or 50%) are statistically indistinguishable between 
the high and low AGN luminosity subsamples, but are 
different at > 5cr between the high and low absorption 
strength subsamples (Table [5]). 

Since we would be dividing on the primary driver behind the 
difference in starburst contribution; the timescale for an outflow 
is thought to be around an order of magnitude shorter than the 
timescale for an AGN, so the outflows would be a pseudorandom 
outcome of a luminous AGN. 
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Fig. 12. — Probability Distribution Function for starburst contribution divided according to their AGN luminosities; in black for the 
whole sample, in red for the objects with Lagn > lO^^'^L©, and in orange for the objects with Lagn < lO^-^'^L©. The difference between 
the tw o PDFs is smaller than in figure [TT] suggesting that the IR luminosity of the AGN is not a good proxy for outfiow strength, see 



Overall therefore, we find that radiatively driven out- 
flows from an AGN with absorption strengths > 3500 km 
s"-*^ act to curtail star formation in their host galaxies. 
We also find that this effect is (at least largely) relative; 
such outflows reduce the contribution from star forma- 
tion to the total IR luminosity to less than ~ 25%. We 
also propose that the infrared luminosity of the AGN is 
not a good proxy for the degree of AGN feedback that 
is taking place. Finally, since the IR properties of our 
sample are consistent with being drawn randomly from 
the reddened QSO population, we conclude that this is 
true generally for reddened QSOs. 

These results are consistent with the idea that star- 
burst and AGN activity crudely correlate with each other 
in active galaxies - a more luminous starburst means we 
are more likely to see a luminous AGl^H - but where 
we add that a radiatively driven outflow can curtail 
the relative luminosity of the starburst on much shorter 
timescales than the typical lifetimes of a starburst or 
AGN. 

We conclude on a cautionary note. Our results appear 
solid, but they are based on a small sample, which makes 

this is consistent with the results i n Table [4] p revious work 
on ULIRGs, see e.g. the luminosities in lFarrah et al.l ((2003 



the errors difficult to estimate. We have used what we be- 
lieve to be a robust error estimation method (jack-knife 
resampling. Tables H] & [5]) , but with only 31 objects we 
cannot reliably measure the error function, since compre- 
hensive resampling methods are not possible. Moreover, 
we cannot test the robustness of the adopted boundary 
for dividing the PDFs of AS2 = 3500 k m s~^ This 
boundary was motivated by the analysis in §4.2.2| but we 
would ideally like to explore the consequences of varying 
this boundary by up to a few thousand km s~^ in either 
direction. We did perform a basic test of this in the ap- 
pendix, by adopting instead a boundary of AS2 = 5000 
km s~^, but we do not consider this to be a robust as- 
sessment of how sensitive our results are to the choice of 
boundary. Given the distribution of AS2 values of our 
sample however, we cannot perform more comprehensive 
tests. A similar argument applies to the choice of lu- 
minosity cut for dividing the AGN PDFs in Figure [121 
and to the choice of template library used to model the 
luminosities (since several alternatives are available, see 
§3.ip . To explore these issues properly would require a 
sample at least a factor of two larger than that used here. 
Aside from sample size, we would also like photometry 
longward of 200/im so as to measure the emission from 
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cold dust heated mainly by star formation, and higher 
quality optical spectra to resolve Mgll kinematics and 
look at the absorption strengths of at least one other 
species. 

5. CONCLUSIONS 

We have presented a sample of 31 FeLoBAL QSOs 
from the SDSS at O.S < z < 1.8. These QSOs have 
broad, deep absorption troughs in their rest-frame UV 
spectra that are unambiguous signatures of radiatively 
driven outflows powered by an AGN. Previous work has 
suggested that FeLoBAL QSOs are IR-luminous, and 
sometimes harbor star formation rates of up to a few 
thousand solar masses per year. Furthermore, there is 
evidence that the AGN-driven outflows can extend up 
to several kiloparsecs into the host galaxies. FeLoBAL 
QSOs are thus ideal laboratories for studying the effects 
of radiatively-driven outflows from an AGN on obscured 
star formation. 

We selected our sample purely on the basis of their 
rest-frame UV spectral properties, and assembled for 
them optical through far-IR photometry from the SDSS, 
2MASS, UKIDSS, WISE and Spitzer. We then fit these 
data with radiative transfer models for the IR emission 
from AGN and starbursts to extract best-fit AGN and 
starburst IR luminosities. We then compared these lu- 
minosities to the strength of their outflows as inferred 
from the absorption properties of the MgIIA2799A line 
in their SDSS spectra. Our conclusions are: 

1 - FeLoBAL QSOs are luminous in the IR. All of our 
sample have total IR luminosities in excess of lO^^L©. 
Nearly one-third of the sample exceed lO^^L©. The dom- 
inant power source behind the IR emission is, in most 
cases, AGN activity. A pure AGN is either the most 
likely power source, or consistent within the 90% confi- 
dence interval, for eleven of the sample. A starburst com- 
ponent is required for the remaining objects, but in only 
twelve of these objects is the starburst more luminous 
than IO^^Lq, and in only three objects is the starburst 
more luminous than the AGN. The mean AGN contri- 
bution to the total IR luminosity is ~ 76%. The spread 
in AGN contribution to the total IR luminosity is how- 
ever wide, spanning 0.2 to 1.0. Overall, the IR properties 
of FeLoBAL QSOs appear consistent with those of the 
general red/dusty QSO population. We do not however 
find convincing evidence that FeLoBAL QSOs are more 
likely to be a ULIRG to QSO transition phase than the 
general red QSO population. 

2 - We find no convincing relationship between the 
strength of the outflows and the total IR luminosity, or 
between the strength of the outflows and either the star- 
burst or AGN luminosities. Conversely, we find a clear 
relationship between the strength of the outflows and the 
contribution from star formation to the total IR luminos- 
ity. If we divide our sample in two at an outflow strength 
boundary of AS2 = 3500km s~^ and construct probabil- 
ity distribution functions for the starburst contribution 
for the two subsamples, then the low absorption strength 
subsample shows a higher chance of a higher starburst 
contribution than the sample as a whole, and a much 
higher chance of a higher starburst contribution than 
the high absorption strength subsample (Figure [TT]). We 
quantify this by extracting the probabilities of obtain- 
ing starburst contributions to the total IR luminosity in 



excess of 25% (Tabled rows 1-3). We find, at > 5<j sig- 
nificance, a higher chance of obtaining a starburst con- 
tribution in excess of 25% in the low absorption strength 
sample compared to the high absorption strength sam- 
ple. 

3 - This anticorrelation between outflow strength and 
the contribution from star formation to the total IR 
luminosity is straightforwardly interpreted as the out- 
flow from the AGN curtailing star formation in the host 
galaxy. There are however several other ways that it 
could arise. The most obvious alternative is that stronger 
outflows reflect an increase in the IR emission from the 
AGN, but have no effect on the starburst; if this is the 
case then the starburst contribution would decline with 
stronger outflows, but without there being a direct re- 
lationship behind the decline. To test this possibility 
we divided our sample into two at an AGN IR lumi- 
nosity boundary of Lagn < lO^^-^L© and constructed 
probability distribution functions for the starburst con- 
tribution for these two subsamples (Figure [T2]). If it is 
the case that outflow strength is a proxy for AGN lu- 
minosity, then we should see a bigger difference between 
these two PDFs than that seen between the PDFs in 
Figure [m Instead, the PDFs show a smaller difference. 
Furthermore, we find that (a) we are (marginally) more 
successful in finding systems with a large starburst con- 
tribution to the total IR emission by selecting on weak 
outflows than we are by selecting on low AGN luminosity, 
and (b) the probabilities of a starburst contribution in 
excess of 25% are statistically indistinguishable between 
the high and low AGN luminosity subsamples, but are 
different at > 5<j between the high and low absorption 
strength subsamples (Table [5]). We considered several 
further alternative possibilities ( §4.2.3p but did not find 
any of them convincing. 

4 - We therefore conclude that strong, radiatively 
driven outflows in FeLoBAL QSOs can have a dramatic, 
negative effect on obscured star formation in their host 
galaxies. This is the most direct evidence yet obtained 
for 'quasar mode' feedback in the QSO population at 
high redshifts. We find that outflows with an absorption 
strength in MgIIA2799A of greater than AS2 = 3500km 
s~^ (Table [3]) act to curtail luminous star formation in 
their host galaxies. We also find that this effect is at 
least largely relative; the starburst luminosity is reduced 
to less than about 25% of the total luminosity. Finally, 
we propose that the magnitude of this effect is not de- 
ducible from the IR luminosity of the AGN. 
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Fig. 13. — Two examples of comparisons between our absorption strength measures, and those for the same objects but by different 
authors. Lejt: AS2 measured by us and by Gibson et al. (2 009). Right: ASq measured by us and bv I Allen et al.l (|201in . In general the 
measures are consistent, albeit with some scatter, see ^4.2.11 
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APPENDIX 

THE EFFECTS OF DIFFERENT ABSORPTION STRENGTH MEASURES 

As described in ^4.2. ![ measuring the properties of BALs is not straightforward. So, we here summarize the checks 
we performed to see if our results are robust against the choice and method of BAL measurement. 

First, we check our measures by con i paring theni to th e three independent measurements already in the literature 
(jTrump et al.ll2QQ6l : iGibson et al .1120 091 : 1 Allen et al.|[2QllL though we note that the different absorption strength mea- 
sures used by these authors together with the evolving reduction of the SDSS spectra means the samples in them 
are not identical). We obtain generally reasonable agreement (Figure [T3|). We conclude that our measurements of 
absorption strength are acceptable. We note though that the scatter in this figure is larger than the formal errors on 
the fits, by a factor of two or more in some cases. This likely reflects systematic uncertainties arising from continuum 
placement. We do not attempt to account for these systematics, or quote them, as they are difficult to estimate 
robustly. Instead, we use our absorption strength measurements only to make comparisons to each other, where they 
should be reliable. As a consequence, we do not quote individual errors on the fits (which are usually of order 100 — 200 
km s~^) as we feel these numbers are misleadingly small. 

Next, we check whether the results in Figure [TOl depend on the adopted parametrization of absorption strength. In 
Figures [14] & [15] we replot these figures with both a stricter and more relaxed definition of absorption strength (ASq 
& AS4, see Table [3]). In all cases we recover comparable results. It is interesting that the relation is slightly offset 
from zero in the bottom right panel of Figure [151 perhaps suggesting that the anticorrelation is not driven by troughs 
much narrower than 2000km s~^. We conclude that our results are reasonably robust to the choice of parametrization 
of absorption strength. 

Finally, we check whether the results in ^4.2.21 depend on the method used to measure the absorption strengths. As 
our sample is small, we cannot robustly test the effects of different techniques to measure absorption strength. We 
can however do two basic checks. First, we keep our AS2 measurements, but change the boundary for the low vs high 
absorption strength samples from 3500 km s~^ to 5000 km s ~^. S econd, we substitute our AS2 measurements with 
the ASo measurements for the same objects in IGibson et al.l ()2009f ). adjusting the boundary to ASq = 4000 km s~^. 
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Fig. 14. — Consistency checks to see if the results in Figure fTOl depend on the parametrization of absorption strength used. The left 
column uses ASq (a strict measure) while the right column uses AS4 (a relaxed measure). Top row: total IR luminosity. Bottom row: 
AGN luminosity. We obtain consistent results in all cases. 

Since [Gibson et ani2QQ9l are more stringent about what constitutes a BAL QSO in comparison to 'Trump et al.| [2QQ6l 
(from where we chose our sample), this is also a test on a smaller, 'golden' sample. We plot the resulting PDFs for 
AGN fractions in Figure [16] and present the confidences on luminosity ranges and AGN fractions in Table [6] for both 
tests. We obtain the same results - an anticorrelation between absorption strength and fractional AGN luminosity - at 
a similar degree of significance. We conclude that our result is robust against the method used to measure absorption 
strengths. 
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Fig. 15. — Further consistency checks to see if the results in Figure fTol depend on the parametrization of absorption strength used. The 
left column uses ASq (a strict measure) while the right column uses AS4 (a relaxed measure). Top row: starburst luminosity. Bottom row: 
AGN fraction. We obtain consistent results in all cases. 
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Fig. 16. — Consistency checks to see if the results in Figure [TT] depend on the parametrization of absorption strength used. The left 
panel uses the same parametrization, but a different bound ary of AS2 = 5 000 k m s"-*^. The right panel substitutes the AS2 measurements 
performed by us with the ASq measurements performed by I Gibson et al.l (|2009l ). We obtain the same results in both cases (see also Table 



